Buckling transition of nematic gels in confined geometry 
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A spontaneous buckling transition in thin layers of monodomain nematic liquid crystalline gel 
was observed by polarized light microscopy. The coupling between the orientational ordering of 
liquid crystalline solvent and the translational ordering of crosslinked polymer backbones inside the 
nematic gel contributes to such buckling transition. As the nematic mesogens become more ordered 
when the gel is cooled down from a higher gelation temperature, the polymeric backbones tend to 
elongate along the direction parallel to the nematic director, which is perpendicular to the rigid 
glass surfaces in the experimental setup. The shape change of such confined gel sample lead to the 
spontaneous buckling of polymeric network and the spatial modulation of nematic liquid crystalline 
director, which is observed as the stripe patterns. The instability analysis was used to explain such 
transitions, and the relationship between the critical field, stripe's wavelength and temperature can 
be explained qualitatively by the rubber elasticity theory for liquid crystalline gels. 
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I. INTRODUCTION 

Liquid crystalline gelsjT] refer to special soft mate- 
rials that incorporate the symmetry properties of liq- 
uid crystalline J2] phases into the crosslinked polymeric 
backbones, thus the translational response of crosslinked 
polymeric networks and the orientational response of liq- 
uid crystalline mesogens are coupled together. Among 
all the possible liquid crystalline gels, nematic gel has 
the simplest symmetry where the crosslinked polymeric 
backbones are spontaneously elongated along one cer- 
tain direction (usually the nematic director n) under 
the effect of symmetry broken properties of the nematic 
solvent. The uniaxial prolate ellipsoidal polymer back- 
bones can be described by a step length tensor = 
+ (in — lx)nirij and the anisotropic parameter r is 
defined as the ratio of effective step length of polymer 
coil parallel (Zy) and perpendicular (ZjJ to the nematic 
director n. The value of r depends on the symmetry 
properties of the nematic solvent : r increases as the sys- 
tem becomes more ordered, i.e. lower temperature for 
thermotropic liquid crystal. When there are no rigid 
mechanical constrains, such relationship can be verified 
by observing the macroscopic shape change of the ne- 
matic gel between the isotropic phase and the nematic 
phase [31 H]. As the temperature decreases, the mono- 
domain nematic gel sample will become elongated along 
the direction of nematic director when sample changes 
from the isotropic phase into the nematic phase. When 
the material is confined by rigid boundaries, for exam- 
ple, on the direction of elongation during the cooling, 
a buckling transition is expected to happen and it has 
been experimentally observed as stripe patterns under 
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polarized light microscopy [5 . Here, we report another 
buckling transition in thin layers of same nematic liq- 
uid crystalline gels within different confined geometries. 
The physical reasons for such buckling transition can be 
qualitatively interpreted by the coupling between the me- 
chanical response of crosslinked polymeric backbones and 
the orientation response of the nematic solvent. The in- 
stability analysis were applied to explain the experimen- 
tal phenomena such as the temperature dependence of 
critical point and wavelength of periodic patterns. The 
study about this buckling phenomena is helpful to pro- 
vide insights to the buckling transition found in other 
soft materials, i.e. microtubules [B] , F-actin networks[7]. 



II. MATERIAL AND EXPERIMENTAL 

Nematic gel material was synthesized in Kornfield's 
group[51 13 EH EQ- Briefly, 5-wt% of ABA triblock 
copolymer, which consist of polystyrene as end blocks 
and side group liquid crystalline polymer as middle 
blocks, were dissolved into a nematic solvent (4-n-pentyl- 
4'-cyanobiphcnyl, 5CB). The formation of weak phys- 
ical network is controlled by the order parameter of 
the solvent: the polystyrene end blocks are soluble in 
the isotropic phase and aggregate in the nematic phase. 
The phase transition temperature of such nematic gel 
(Tni ~ 37°C) is very close to the transition tempera- 
ture of 5CB (Tni ~ 35°C), and the reversibility of the 
physical crosslinking mechanism allows repeatable exper- 
iments being easily conducted on the same sample. 

The nematic gel is loaded in a 25/im thick homeotropic 
electro-opitcal cell, in which a thin layer of ra,n-dimethyl- 
n-octadecyl-3-aminopropyl-trimethoxysilyl-chloride were 
spin coated on the surface of transparent indium-tin ox- 
ide conductors of glass slides. In the presence of a strong 
applied electric field (E = 3V/fim, 1kHz) across the 
cell, the nematic mesogens can be easily aligned verti- 



FIG. 2: Experimental measurement about the temperature 
dependence of | (a) | wavelength of stripes in the buckled state 
and (b) critical field (Ec) of the buckling transition in nematic 
gel. 



FIG. 1: Optical micrographs of birefringent stripe pattern of 
buckled nematic gel observed through polarized optical mi- 
croscopy with applied electric field decreased to zero at dif- 
ferent temperatures. The scaling bars stand for 50/im in all 
the images. 



cally throughout the cell, where the long axis pointing 
perpendicularly to the boundary surfaces. The temper- 
ature of the sample was controlled by a peltier-based 
microscope stage during the observation. Initially, the 
sample was heated up to 45°C in the isotropic phase with 
no crosslinked polymer networks. The sample was cooled 
down (2°C / minute) across its Tni to certain final temper- 
atures (Tf) and a thin layer film of mono-domain nematic 
gel was obtained during the gelation throughout the cell 
volume, in which n points perpendicularly to the bound- 
ary surfaces. The sample appeared homogenous dark 
under the cross polarized optical microscopy while the 
aligning field maintained its original magnitude. When 
the electric filed was turned off, birefringent stripe pat- 
terns with wavelength about 5/xm appeared throughout 
the sample, as shown in Fig. [T] Both the wavelength of 
stripe pattern and the critical field (Ec), at which the 
sample changes from homogeneous dark to birefringent 
patterned, depend on the sample's final cooling temper- 
ature (Tf), such temperature dependance are recorded 
and plotted in Fig. [2] It can be seen that both Ec and 
wavelength stay in a plateau when 10°C< Tf <24°C , and 
when Tf > 24°C Eq decreases as Tf increases, while the 
wavelength increases. 



III. DISCUSSIONS 

The driving force of such transition can be attributed 
to the thermo-mechanical-optical coupling between 
the nematic liquid crystalline solvent and anisotropic 
crosslinked polymeric backbones in the nematic gel 
within a confined boundary condition. Diagrams in Fig. [3] 
are used to interpret the physical reasons for such tran- 
sition. Initially, the monodomain nematic gel is formed 
at a higher initial temperature (T Tni) within the 
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FIG. 3: Diagrams of the buckling transition in nematic gel. 
The rigid boundaries cause the material buckling within it- 
self as the polymeric backbones elongate in a more ordered 
state |(b)| in comparison to a less ordered gelation state |(a)| 
at which the liquid crystal mesogens are aligned vertically by 
the applied electric field E. 



glass cell as the e lectric field is applied across the cell, as 
show in Fig. |3(a) The micro picture of the nematic gel is 
sketched as the ellipsoid and macroscopic shape of the gel 
is expressed as the square shape in the diagram. As the 
temperature is lowered to final temperature (Tf < T), 
the system becomes more ordered and the anisotropic 
polymeric coil will elongate along the nematic director's 
direction, as shown in the ellipsoid of Fig. |3(b)| If there 
is no boundaries to confine the shape of the material, 
the macroscopic shape of the mono-domain gel sample 
will elongate vertically, which is illustrated as the rect- 
angular dashed line. Such "artificial muscle" effect has 
been observed experimentally in nematic elastomers and 
gelsjl]. When the sample is put into an environment with 
rigid constrains, e.g. a cell with two glass slides glued to- 
gether, the material has to buckle within the boundaries 
in order to gain the elongation along the direction per- 
pendicular to the rigid boundaries. Due to the coupling 
between the translational response of the polymeric coil 
and rotational response of the liquid crystalline mesogen, 
the nematic director n will rotate correspondingly. The 
applied electric field with enough magnitude can be used 
to keep the nematic director aligned vertically, and the 
material will buckle spontaneously as the electric field is 
decreased due to the instability behavior. To prove our 
physical explanation, detailed analytical calculations are 
conducted and discussed in the following. 
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A. Modeling and Free Energy Calculation 

The coordinate origin is selected at the middle of the 
cell gap and the z-axis is perpendicular to the cell bound- 
aries (z = ±d/2). Initially, the aligning electric field 
is applied along the z-axis (E = Eqz) and the nematic 
director is initially aligned vertically (n° = z) as well, 
as shown in Fig. 3(a)| The superscript is added onto 
parameters for the initial gelation state of the material. 
The polymeric networks are formed during the gelation 
process at temperature Tj, with r° as the anisotropic pa- 
rameter. When the material is cooled to a lower final 
temperature Tf < Tj, the crosslinked polymeric network 
will be more elongated along nematic director n as the 
nematic solvent become more ordered and the anisotropy 
parameter r that is larger than the initial value, r > r°. 
When the electric field is decreased to zero (E = 0), the 



elongation of the polymeric coils within the nematic liq- 
uid crystalline gel buckle with a displacement field within 
the xz- plane: R = Q cos kx cos gzx + rjsmkxsvaqzi, in 
which k is the wavevector of the shear wave within xy- 
plane, and we select it along the x-axis, q is the wavevec- 
tor along the z-direction. The value of q is determined 
by the sample thickness d as q = ir/d for the first har- 
monic mode. C an d f] are related by the incompress- 
ibility conditionQQ as qr\ = fc£. Such shear motion of the 
polymeric network will induce the nematic director to ro- 
tate within xz-plane of small amplitude £ about z-axis: 
n = £ cos kx sin gzx + z. 

We plug these conditions into the formulas for free en- 
ergy density, which include the Frank curvature elastic 
energy of the nematic solvent [5] and the nematic rubber 
elastic energ of crosslinked polymeric backbones pQ: 



F = - (K s (V • n) 2 + K B (n x (V x n)) 2 ) - -e a (E • n) 2 

+ \v ■ Tr(t • A T • F 1 • 4) - l -ApL • Tr(4 T • n • A_ - nP • A_ T • n ■ A_). (1) 



In Eq. [T] Ks and Kb are the curvature elastic constants 
for the nematic solvent; fi is the shear modulus of the 
gel and A is the semisoftness coefficient; A is the Cauchy 
strain tensor, Ay = 5y + dRjdxj. 



It can be seen that the free energy is proportional to the 
square of the perturbation's amplitude £ 2 . If / < 0, the 
perturbed state would be stable and the system would 
have a transition with £ ^ 0; on the other hand, if / > 0, 
the perturbed state is unstable and the system would stay 
in its initial state with C, = 0. / = is the critical point 
at which the transition starts. In this way, the values 



I 

By averaging this free energy density over the space, 
minimizing with respect to £, and keeping only terms 
only to second order in £, the free energy density / in 
the material can be written in Eq. [2] 



(2) 
I 

of the material's physical parameters determine the in- 
stability behavior under different external experimental 
conditions, e.g. applied field, temperature. We can plug 
the known parameters of the nematic liquid crystalline 
gelOOEEI] into Eq.|2] as d = 25/xm, fi = 50Jm- 3 , A = 0.1, 
K B = lQ- n 3m-\K s = 1.5 x lO^Jm- 1 , e a = 15e . 
We choose r° = 1.5 for the initial gelation state at tern- 
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/<! x 

Vr(l - rO + r(A - 1 + ^-+K s ^+K Bq ^ + r0 ^ 

( k 6 r(l + Ar) K ^+ q 4 r°(r + (A - 1 + ^ + K B<? y2 + , _ - on 
V A* M 
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perature Tj, and r depends on the order parameter of by setting k = 0: 
the nematic solvent in the gel for specific temperature 
Tf with elongated polymeric coil. 



B. Instability Analysis: Critical Field 

For the case of a homogeneous birefringence change as 
the basic mode of transition, / can be further simplified 



220 / r + (r - l)r° + r 2 ( A - 1 + 



fj,( q r 
4r 



ft 



r A-l 



e a E 2 +K B q 2 



(3) 



Fig. |4(a) shows the plot of /fc_>o as a function of electric 
field's intensity E, which corresponds to the situation 
when we decrease the electric field at a certain temper- 
ature of final state. The free energy is positive when 
the electric field is still large enough to keep the vertical 
alignment of the nematic directors across the cell, the 
buckling transition is energetic unfavorable; as the elec- 
tric field is further decreased, the energy become nega- 
tive, which means the system transforms into the buckled 
state to minimize the free energy. The critical electric 
field Eq can be found at the point where the free energy 
equals zero. The relationship between the Eq and r can 
be studied numerically and plotted in Fig. |4(b)[ where 
Ec increases with r, in another word, Eq decreases with 
final temperature TV. This agree qualitatively with our 



experimental measurement in Fig. |2(b)| when r = 3.5, 
critical field is calculated as Eq = 0.48V//im comparing 
with experimental value 0.56V//im at 31°C; when r = 5.8 
corresponding to a lower Tf, critical field is calculated 
as Eq = 0.621V//Ltm comparing with experimental value 
0.852V//zm at 25°C. 



C. Instability Analysis: Stripe Wavelength 

Since the nematic gel buckles when the applied electric 
field removed (E = 0), the free energy / depends on the 
wavevector k in xy-plane: 



Aq 2 r{\ - r° + r(A - 1 + 

t 



K s k 2 + K B q 2 
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+k\ (1 + r°- Ar° +^(K B + K s r )) + k 2 q 2 ( (3 - r°)r° + r(l + r°(3A - 6 + + r )) 
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Fig. 4(c) shows the plot of Je^q as a function of wavevec- 
tore k. It can be seen that free energy is negative at 
a finite k, which corresponding to the stripe pattern 
observed experimetally. The spatial modulation of the 
translational order (R) and the orientational order (n) 
further minimize the free energy, which is similar to the 
stripe pattern observed in previous experimental observa- 



tion in the planar aligned sample 5 . Furthermore, wave- 
length can be numerically calculated as a function of the 
sample's anisotropic parameter r at final state, which 
is plotted in Fig. |4(d)| It can be seen that wavelength 
decreases as r increases, in another word, stripe's wave- 
length is smaller for lower final temperature Tf. This 
agrees with our experimental observations in Fig. 2(a) 
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length with temperature. The theoretical calculation can 
only be compared with experimental data qualitatively. 
It can be seen that it agrees well qualitatively with the 
experimental measurements. 




IV. CONCLUSIONS 



(a)/fc_ vs - E 



(b)Ec vs - r 




(c)/b—o vs - k 



(d)A vs. r 



FIG. 4: The fr ee e nergy density (/) depends both on (a) 
electric field and (c) wavevector k. The free energy is mini- 



mized by decreasing applied field to zero and spatial modu- 
lation of director field n with a finite k. Both critical field 
Eq and wavelength depend on the anisotropic properties of 
final buckled state: as the nematic gel is more ordered a t its 
final state (corresponding lower Tf), (b) Ec increases and (d) 
wavelength decreases. 



when r — 3.5, the wavelength is calculated as 6.48/im 
comparing with experimental value 6.46/im at 31°C; when 
r = 5.8 corresponding to a lower Tf, the wavelength is 
calculated as 3.82/mi comparing with experimental value 
3.82^im at 25°C. 

Currently, the analytic relationship between tempera- 
ture T (or the order parameter of nematic liquid crystals) 
and the anisotropy parameter r of polymeric networks is 
not known experimentally. Therefore we can not fit our 
experimental measurement of the critical field and wave- 



In summary, the spontaneous buckling transitions 
of thin layers of nematic liquid crystalline gel in a 
homeotropic cell were observed by polarized light mi- 
croscopy. This is good example to show the coupling be- 
tween the liquid crystalline ordering and the crosslinkcd 
polymer backbones inside the nematic gel material. As 
the nematic mesogens become more ordered when the gel 
is cooled down from the initial crosslinking stage with a 
higher temperature, the polymer network tends to elon- 
gate along the direction parallel to the initial nematic 
director, which is perpendicular to the rigid glass sur- 
faces in the experimental setup. The shape change of 
such confined gel sample lead to the spontaneous buck- 
ling transition. The applied electric field will change the 
in stability behavior, and spatial modulated stripe pat- 
tern in orientational ordering of nematic solvent helps 
to accommodate the buckling transformation of gel net- 
work and minimize the free energy. The experimental 
observation and measurement can be can be explained 
qualitatively at different temperature. 
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